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Electronic structure calculations within screened-exchange hybrid density functional theory show 
that Cui n Ga antisites in both CuInSe2 and CuGaSe2 are localized hole traps, which can be attributed 
to the experimentally observed N2 level. In contrast, Gac u antisites and their defect complexes 
with copper vacancies exhibit an electron trap level, which can limit the open-circuit voltage and 
efficiency in Ga-rich Cu(In,Ga)Se2 alloys. Low-temperature photoluminescence measurements in 
CuGaSe2 thin-film solar cells show a free-to-bound transition at an energy of 1.48 eV, in very good 
agreement with the calculated transition energy for the Gac u antisite. Since the intrinsic DX center 
InDx does not exhibit a pinning level within the band gap of CuInSe2, metastable DX behaviour 
can only be expected for Gac u antisites. 



Thin-film solar cells with aborbers based on 
Cu(In,Ga)Se2 alloys currenly achieve record efficiencies 
of 20.3% and represent one of the most promising 
technology for large-scale industrial production pQ. 
Efficiencies up to 14.2% can also be obained with the 
ternary boundary phase CuInSe2, whereas the efficien- 
cies of pure CuGaSe2 cells are still limited to below 10 
% 0. Intrinsic point defects acting as recombination 
centers are likely to limit the open-circuit voltage and 
therefore the solar cell efficiency of these absorbers. 
The point defect physics of these chalcopyrites has 
been extensively studied by experimental methods, 
such as electrical and optical spectroscopy methods, 
and theoretical approaches, mostly calculations based 
on density functional theory [3Tll3j. A hole trap level 
in the range between 0.15-0.35 eV, often named N2, 
has been observed using admittance and deep-level 
transient spectroscopy (DLTS) and Hall measurements 
in CuInSe 2 [HHS], CuGaSe 2 and Cu(In,Ga)Se 2 

[14) , H"5 ] [20H22] b A measured activation energy between 
0.05-0.20 eV, has been attributed to an interface defect 
in Cu(In,Ga)Se2 and was named Nl [20 , although this 
denomination is ambiguously used and it is unclear 
wether the Nl response is associated to a defect at 
all [23]. Indeed, various metastable effects have been 
oberserved in CIGSe devices such as persistent pho- 
toconductivity [23], the increase of the open-circuit 
voltage upon white-light soaking [25 , an increase of the 
space-charge upon illumination [26 or reverse-biasing 
[27] accompanied with a decrease of the fill factor [28] as 
well as capacitance relaxation on long time scales after 
light-soaking [25] . 

Based on electronic structure calculations within local 
density functional theory, two intrinsic point defects and 
their complexes with copper vacancies have been pro- 
posed to exhibit metastable properties in CuInSe 2 and 
CuGaSe2: the intrinsic indium and gallium DX centers 
(In, Ga)ox[in] [30] and the selenium vacancy Vs e [8 J or se- 
lenium vacancy-copper vacancy complex Vs e — Vc u [S] • 



These theoretical results have been invoked [3TH51] to ex- 
plain the experimentally observed light and voltage-bias 
induced metastabilities. However, metastable point de- 
fects are not the only possible explanation. Copper mi- 
gration in the space charge region [551 13"6"] . deep acceptor 
levels in the CdS buffer layer[37], an electron-injection 
dependent barrier at the molybdenum back-contact of 
the device [23] . or the presence of a p+- layer in conjunc- 
tion with a shallow donor level at the buffer absorber 
interface [38] have all been put forward as possible ex- 
planations. Therefore, metastabilities in Cu(In,Ga)Se 2 
based devices and their possible relation to the N2 and 
Nl levels remain puzzling and it is not clear whether a 
single explanation is sufficient to explain all of the ob- 
served phenomena [23] [39]. Despite extensive efforts, no 
concise picture of the point defect physics that matches 
with all of the experimental findings has yet emerged. 

In this letter, we show that Cuj n Ga antisites are local- 
ized hole traps in both CuInSe 2 and CuGaSe 2 , which can 
be attributed to the experimentally observed N2 level, 
while Gacu antisites and their defect complexes with 
copper vacancies exhibit an electron trap level, which 
can limit the open-circuit voltage and efficiency in Ga- 
rich Cu(In,Ga)Se 2 alloys. Our results are based on elec- 
tronic structure calculations within screened-exchange 
hybrid density functional theory and are put in context 
to low-temperature photoluminescence measurements of 
CuGaSe2 thin-film solar cells. 

We have carried out hybrid density functional calcula- 
tions using the HSE06 functional [301 HI] as implemented 
in VASP[42] with an adapted exchange-screening param- 
eter of 0.13 A _1 [43], which allows to closely match the 
experimentally observed band gaps for both CuInSe 2 , 
CuGaSe2 and other chalcopyrite phases [IT]. The ap- 
proach allows to overcome the band-gap problem, to di- 
rectly analyze defect levels within the gap and improves 
the description of localized and correlated copper d elec- 
trons O [45]. Thus, using the HSE06 functional will 
lead to more accurate formation enthalpies of point de- 
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FIG. 1. Stability diagram for CuInSe2. The defect formation 
enthalpies in Fig. [2] are discussed in terms of the chemical po- 
tential at points A,B and C. The stability range for CuGaSe2 
(not shown) has the same shape, but a somewhat larger ex- 
tent due to its higher formation enthalpy (CuInSe2: -2.37 eV, 
CuGaSe 2 : -2.67 eV). 



fects than LDA calculations. A 2x2x2 r-centered k-point 
grid has been used for supercells with 64 atoms as well as 
with 216 atoms. Ion positions were relaxed until forces 
were converg ed to below 0.05 and 0.1 eV/A, for super- 
cells of 64 and 216 atoms, respectively. The calculation 
of supercells of 216 atoms with a 2x2x2 k-point grid are 
computationally extremely costly, but necessary in or- 
der to observe unambiguously localized defect levels of 
Cuj n Ga and Ga Cu . All reference phases presented in the 
stability diagram (Fig. [IJ were calculated using the same 
functional. The point defect formation enthalpies were 
calculated as function of the chemical potentials of the 
constituents A/x, referenced to the elemental phases and 
the Fermi energy ep according to the common formula 
as e.g. in Ref. [7 . The potential alignment and im- 
age charge corrections have been carefully carried out as 
described in Ref. |46l 

Since high-efficiency Cu(In,Ga)Se2 absorber material 
is prepared under a highly selenium-rich atmosphere, 
e.g. with a selenium to metal flux ratio of 5 [47], 
it is instructive to interpret the defect physics for 
this material under selenium-rich conditions on the Se- 
Cu(In,Ga)Se2 phase boundary in the calculated stabil- 
ity diagram (point A in Fig. [IJ. In contrast to a previ- 
ously calculated stability diagram [7], our calculations 
show a phase boundary between Cu2Se and CuInSe2, 
which is in line with the experimental phase diagram 
and the observation of Cu2Se precipitates under certain 
processing conditions [48 j (point B). For Cu-poor com- 
positions, which yield the highest conversion efficiency, 
Cu(In,Ga)Se2 is a highly compensated semiconductor. In 
this case, the charge neutrality condition and thus the 
Fermi-energy is essentially determined by the concentra- 
tion of donors and acceptors with the lowest formation 
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FIG. 2. Defect formation enthalpies, charge transition lev- 
els and the determined intrinsic Fermi levels in CuInSe2 and 
CuGaSe2 for chemical potentials corresponding to points A,B 
and C in Fig. [Tj 



energies (see Fig. [2|. For the chemical potentials at 
points A, B and C the material turns out to be p-type, 
while it becomes n-type for maximal Cu- and In-rich con- 
ditions (not shown) . Fig. [2] shows the calculated defect 
formation enthalpies for the various chemical potentials. 

One of the most intriguing result is that Cui n .Ga &n- 
tisites in both materials, CuInSe2 and CuGaSe2, can 
have equally low formation enthalpies as copper vacan- 
cies and thus also act as compensating defects. This 
finding is consistent with large concentrations of Cui n in 
CuInSe2 recently reported using wavelength dispersive x- 
ray diffraction even for copper-poor material [49 . How- 
ever, when the chemical potentials are shifted towards 
metal-rich conditions (e.g. point C), it is seen that the 
formation enthalpy of Cui nj Ga at the intrinsic Fermi level 
(vertical arrow) increases, while it does not change much 



for (In, Ga) Cu and V, 



'Cu- 



In° 



The analysis of the density of states of Cuj n Ga and 
reveals an empty narrow defect band above the 



Cu.Ga 



VBM for Cui n in CuInSe 2 (at 0.27 eV) and Cu Ga in 
CuGaSe2 (at 0.32 eV), which can trap two holes (the hole 
density of the empty single-particle defect state of Cuc a 
as obtained from the calculation is displayed in Fig. |4| . 
Furthermore, a localized electron trap level emerges for 
Ga Cu at a single-particle energy of 1.17 eV above the 
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FIG. 3. Aligned density of states of CuGa , Ga^ u and 
(GaQ u — 2Vcu)~ 2 in CuGaSe2 in supercells of 64 and 216 
atoms as obtained with the adapted HSE06 functional. 

VBM in CuGaSe2 (see also Fig. [1]). Forming defect com- 
plexes such as Gacu — Vq u and (Gac u — 2Vc u ) -2 , does 
not affect the position of the defect level. The defect 
level of In Cu is found to be resonant within the CB for 
CuInSe 2 at 1.48 eV and for CuGaSe 2 at 1.46 eV approx- 
imately independent of gallium content. All individual 
single-particle defect levels are approximately constant 
on an absolute energy scale, when comparing their posi- 
tion in CuInSe2 and CuGaSe2 relative to the respective 
VBM. However, the defect levels of In Cu and Ga Cu as 
individual defects do not align on an absolute scale, the 
In Cu defect level being 0.29-0.41 eV higher than the one 
of Ga Cu in both CuInSe2 and CuGaSe 2 as host material. 

Here it should be noted that defect calculations on 
Cu(In,Ga)Se2 have so far been performed with supercells 
smaller or equal to 64 atoms. However, very disperse 
defect bands appear within the gap for the cells with 
64 atoms (see Fig. [3]), which indicates significant self- 
overlap of the defect wavefunctions. If local or semilocal 
density functionals are used, derealization also occurs in 
large supercells up to 216 atom, which is why localized 
deep hole traps have not been identified in the past. This 
proves that an accurate nonlocal treatment of exchange 
and correlation and large supercells are crucial for ob- 
taining the correct localization behaviour of the Cuj n Ga 
defect. 

The Cui n .Ga antisites localize holes, are abundant un- 
der typical preparation conditions (up to 10 20 cm~ 3 at 
850 K deposition temperature and the calculated ther- 
mal transition energy for the process CuQ a +hy B — > Cu Ga 
of 0.20 eV agrees with experimental measurements (e.g. 
0.1-0.3 eV in Ref. (22] and 0.1-0.2 eV in Ref. [21]). Given 
these evidences it is safe to conclude that the N2 hole 
trap level is due to the Cui n .Ga antisite. The fact that 
this level does not occur in all samples can be explained 




FIG. 4. Defect level electron charge density of GaQ u (left: iso- 
surface 0.03 eA -3 ) and hole density of Cug a (right: isosurface 
0.02 eA -3 ) as obtained in supercells with 216 atoms. 



with differing formation enthalpies relevant for different 
preparation conditions (compare points A,B,C in Fig.[2|. 

The Ga Cu antisite, in contrast, shows a clearly local- 
ized electron trap level at a single-particle energy of 1.17 
eV above the VB in CuGaSe 2 and at 1.07 eV in CuInSc 2 
very close to the CBM. Therefore this antisite defect be- 
comes increasingly deep when Ga is alloyed into CuInSc2 , 
due to the rising CB. Since Ga Cu is expected to occur in 
large quantities in Cu(In, Ga)Se2 due to its low formation 
enthalpy, it may limit solar cell efficiency when Ga alloy- 
ing is used to increase the band-gap: When the CBM 
is raised above the position of the Gac u defect level, 
recombination through this defect may limit the open- 
circuit voltage. Since the associates of Gac u with copper 
vacancies display the same single-particle defect level as 
the non-complexed antisite, the complexes may cause the 
same limitations. With regard to the defect complexes 
(In, Ga)c u — 2V Cu we find that their formation enthalpies 
are not particularly low (Fig. [2]). In fact, the total bind- 
ing energy of the complexes with respect to the isolated 
charged species is only about 0.3 eV. 

It might be tempting to relate the experimentally ob- 
served Nl signature to the normal configuration of the 
Gacu defect in Cu(In, Ga)Se 2 . However, we note that 
such an assignment is contradicted by the observation of 
Nl in pure CuInSe2 and the fact that its activation en- 
ergy does not change with increasing Ga-content [5T] [52] . 

To confirm our theoretical results for the Gac u de- 
fect, we conducted temperature-dependent photolumi- 
nescence measurements on a CuGaSe 2 thin-film solar 
cells prepared by a three-stage coevaporation process, as 
used for high-efficiency chalcopyrite solar cell devices [2]. 
For the thin film absorber investigated, a ratio [Cu]/[Ga] 
= 0.87 was measured by x-ray fluorescence analysis 
and the accompanying solar cell showed a device effi- 
ciency of 7%. Photolumincscence (PL) was measured 
using a 670 nm diode-laser as excitation source and 
a thermoelectrically-cooled InGaAs array coupled to a 
0.5 m spectrograph for luminescence detection, with the 
sample placed in a closed-cycle helium cryostat. A pho- 
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FIG. 5. Photoluminescence spectrum at T=15 K measured 
for CuGaSe2 thin films. The black solid line indicates a fit 
consisting of three gaussian-shaped transitions at El = 1.6 
eV, E2 = 1.48 eV and E3 = 1.17 eV (dashed lines). 



toluminescence spectrum obtained at 15 K is shown in 
Fig. [5j Three recombination peaks located at 1.6 eV, 
1.48 eV and 1.17 eV can be clearly distinguished. The 
temperature and excitation intensity dependence of these 
peaks is consistent with an assignment of the peak at 1.6 
eV to a tail-to-band transition, commonly observed in 
Cu-poor chalcopyrites [S3], and the assignment of the 
peak at 1.48 eV to a free-to-bound transition. The lat- 
ter transition energy is in excellent agreement with the 
calculated optical emission energy for the recombination 
process GaJ u + h^ B -> GaJ+ in CuGaSe 2 (1.44 eV) [54] . 
The temperature and excitation dependence of the third 
broad peak observed at about 1.17 eV is consistent with 
a donor-acceptor transition. A possible candidate for this 



transition is the process GaJ u + Cu° 



Ga 



Cu, 



Ga' 



i.e. the recombination of a single electron localized on 
a Gacu antisite with a neighbouring Cuc a hole trap by 
radiative tunneling, which has a calculated transition en- 
ergy of 1.02 eV. 

For metastabilities originating from intrinsic DX cen- 
ters to occur, it is necessary that a DX pinning level 
exists within the gap|10j. From our results (Fig. [2|, 
we conclude that such a pinning level only occurs for 
Ga Cu antisites in CuGaSe 2 (^gx^ 2 = 1.16 eV), but 
not for In Cu antisites in CuInSe 2 (E^^ 3 = 1.31 eV, 
well above the CB). Therefore, Fermi-level pinning and 
metastable effects due to intrinsic DX centers may occur 
in larger band gap Cu(In,Ga)Se2 materials, but not for 
the ternary CuInSe 2 . In order to assess the energy dif- 
ferences responsible for the different DX pinning levels 
as compared to Ref. |10) . it is instructive to compare 
the uncorrected formation enthalpies of In^J and Innx in 



CuInSe 2 within our approach to the ones obtained from 
LDA applying only the static +U valence band correction 
of Ref. |10j. which we were able to reproduce. The uncor- 
rected formation enthalpy of In^ within our approach is 
0.4 eV lower, while the one of Innx is 0.32 eV higher. 
These energy differences, which result from the different 
treatments of exchange and correlation within HSE06 as 
compared to the LDA, which result mostly from the im- 
proved description of the Cu d electrons, directly cause 
the change in the DX pinning levels. 

We have also investigated Vs c — Vc u vacancy pairs, 
which have been held responsible for a variety of metasta- 
bility phenomena in Cu(In,Ga)Se 2 . With respect to this 
defect pair our hybrid functional calculations yield com- 
parable charge transition levels (+/-) and metastable 
relaxation behavior as previously found in LDA-based 
calculations [9 . However, the formation enthalpies of 
this defect complex are higher than 2 eV at the relevant 
Fermi levels at points A,B and C in Fig. |2j). Thus, un- 
der thermal equilibrium conditions, Vs c — Vc u associates 
should occur only in minor quantities (below 10 12 cm -3 
at 850 K deposition temperature as estimated from the 
formation enthalpies). Thus, these defects can only cause 
metastable phenomena if the material is prepared under 
far-from-equilibrium conditions. 

In conclusion, we have shown that Cui n ,Ga defects 
create a localized hole trap level in both CuInSe 2 and 
CuGaSe 2 . Given the large amount of experimental evi- 
dence for a hole trap level in the range 0.15-0.35 eV, we 
conclude that Cui n ,Ga defects should be assigned to the 
N2 level. Ga" u and its complexes display a localized de- 
fect level within the gap in CuGaSe 2 and CuInSe 2 . These 
defects are thus a likely cause for the limited efficiency of 
CuGaSe 2 based wide band devices. Since In^ x does not 
exhibit a Fermi pinning level in CuInSe 2 , metastable DX 
behaviour can only be expected for Gac u antisites. 
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